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Abstract 

Complimenting recent work on the effective field theory of cosmological large scale 
structures, here we present detailed approximate analytical results and further pedagogical 
understanding of the method. We start from the collisionless Boltzmann equation and 
integrate out short modes of a dark matter/dark energy dominated universe (ACDM) 
whose matter is comprised of massive particles as used in cosmological simulations. This 
establishes a long distance effective fluid, valid for length scales larger than the non-linear 
scale ~ lOMpc, and provides the complete description of large scale structure formation. 
Extracting the time dependence, we derive recursion relations that encode the perturbative 
solution. This is exact for the matter dominated era and quite accurate in ACDM also. The 
effective fluid is characterized by physical parameters, including sound speed and viscosity. 
These two fluid parameters play a degenerate role with each other and lead to a relative 
correction from standard perturbation theory of the form ~ 10~^c^ k'^ /H'^ . Starting from 
the linear theory, we calculate corrections to cosmological observables, such as the baryon- 
acoustic-oscillation peak, which we compute semi-analytically at one-loop order. Due to the 
non-zero fluid parameters, the predictions of the effective field theory agree with observation 
much more accurately than standard perturbation theory and we explain why. We also 
discuss corrections from treating dark matter as interacting or wave-like and other issues. 
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1 Introduction 



An effective field theory is a description of a system that captures all the relevant degrees of 
freedom and describes all the relevant physics at a macroscopic scale of interest. The short 
distance ("UV") physics is integrated out and affects the effective field theory only through 
various couplings in a perturbative expansion in the ratio of microphysical UV scale/s to the 
macroscopic ( "IR" ) scale being probed. This technique has been systematically used in particle 
physics and condensed matter physics for many years (e.g., see [H El E]), but has not been 
fully used in astrophysics and cosmology. The large scale properties of the universe acts as an 
important application and is in need of careful analysis. 

Of current fundamental importance is to understand the initial conditions, contents, evolution 
and formation of the universe. It appears to be adequately described by the so-called ACDM 
cosmological model in which the matter content of the universe is primarily dark matter and the 
late time dark energy is adequately described by a cosmological constant. The early universe 
was dominated by a cosmic plasma in which baryons were tightly coupled to photons leading 
to so-called baryon-acoustic-oscillations. The evidence for this model comes from a range of 
sources, including CMB data, lyman-a forrest, curvature constraints, supernovae type lA, weak 
lensing, and (of particular importance to the current discussion) structure formation, etc. In 
this cosmological model, structure formation is primarily driven by the gravitational attraction 
of dark matter, which led to the clumping of baryons including stars, galaxies, and clusters of 
galaxies. This arose from gravitational instabilities of the initial linear density fluctuations that 
were approximately adiabatic, scale- invariant, and Gaussian (e.g., see [H 13 El [7]). 

The power spectrum of large scale structure at late times is corrected from the initial linear 
input in interesting and important ways. For instance, non-linear effects alter the shape of the 
baryon-acoustic-oscillations in the power spectrum. The baryon-acoustic-oscillations are a vitally 
important probe of dark energy as they provide a standard ruler to constrain the cosmological 
expansion history (e.g., see [H |9l [TOl [H]). In general, one needs a proper understanding of the 
departures from linear theory in order to constrain fundamental physics, such as dark energy, 
primordial non-Gaussianity, and other probes of microscopic physics. 

There has been a substantial amount of work in the literature to understand non-linear 
structure formation in the form of cosmological perturbation theory, including [121 [ISl [El [IS 
[l6l[I71[T8l[l9l[20l|2ll[2a[23l|2l|25l[2^ This includes what is usually referred to 

as "standard-perturbation-theory" (SPT). In this approach, the continuity and Euler equations 
for a pressure-less and non-viscous dark matter fluid (with vanishing stress-tensor) is assumed. 
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These non-linear equations for the dark matter are solved perturbatively around the linear 
solution and corrections are obtained order by order, usually truncated at the one-loop order. 
The theory involves integrating /c- modes in the entire domain < /c < oo, which involves treating 
all /c-modes as perturbative. 

The analysis of the current paper stems from the fact that this "standard" procedure neces- 
sarily has a qualitative and quantitative problem. The density fluctuations are not perturbative 
beyond a scale /cnl, the "non-linear scale"; the scale at which density fluctuations are 6 ~ 0(1), 
where gravitational collapse may occur. This scale is roughly Xnl ^ lOMpc or so. Hence there 
are two regimes: k < /cnl which is weakly coupled and perturbative, and k > /e^l which is 
strongly coupled and non-perturbative. For cosmological purposes, we are normally interested 
in the low /^-regime. However, these 2 regimes are coupled by non-linearities, so we must be 
very careful in attempts to describe the low fc-regime perturbatively. The rigorous and complete 
method to do this is that of effective field theory. The procedure is to introduce some arbitrary 
cutoff A on the /c-modes of the fluid. We take this cutoff to be A < /cnl? so that all modes of the 
fluid are perturbative. This means that the high /c-modes {k > A) must be "integrated out". 
In practice this means that these UV modes generate higher order derivative and non-linear 
corrections to the fluid equations for the low /c- modes {k < A). We show that this includes 
terms such as pressure and viscosity; precisely the terms that are assumed to vanish in SPT. 
These terms are a real property of the dark matter fluid and they alter the power spectrum in 
an important and measurable way. These fluid parameters can be determined by matching to 
the full UV theory, i.e., N-body simulations. This furnishes an effective field theory for dark 
matter. This is a fluid that only involves weakly coupled modes and is connected to the full 
microphysical theory through these higher order operators. 

Important earlier work on this topic was performed in the very interesting Ref. [31], where 
this basic conceptual foundation was laid out with particular focus on the issue of back-reaction 
at the scale of the horizon. In our recent work, which we continue here, we (i) focus on sub- 
horizon scales, (ii) obtain an explicit measurement of fluid parameters, (iii) perform an explicit 
computation of the power spectrum, and (iv) provide further insight and clarifications. In the 
present paper we compliment and extend our important recent work in Ref. |32] in which the 
measurement of fluid parameters was performed and the basic framework was put together. 
Here we develop and describe in detail this effective fluid description of dark matter (and by 
extension; all matter, since baryons trace dark matter on large scales), which is the complete 
description of large scale structure formation. In particular, we recapitulate how to extract 
various fluid parameters from N-body simulations and then solve the effective fluid theory to 
some desired order in a perturbative expansion which we formulate recursively. We show how 
to approximately extract the time dependence in a ACDM universe, which connects in a simple 
and intuitive way with the previous standard perturbation theory, but highlights the essential 
differences arising from the fluid parameters. This leads to convenient and quite accurate results. 
Our basic method and key results are summarized in the following discussion: 

If we assume that the matter content of the universe is dominated by non-relativistic matter. 
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primarily dark matter evolving under Newtonian gravity, we can smooth the corresponding 
coUisionless Boltzmann equation in an expanding FRW background. This generates the usual 
continuity and Euler equations. An important point stressed in Refs. [3T| 132] is that the latter 
includes an effective stress-tensor [t^-']k that is sourced by the short-modes 5s. By defining the 
effective stress-tensor by its correlation functions, it can be expanded in terms of the long density 
5i and velocity v\ fields as 



Ha 



+ . 



The individual parameters and cf, = c^^ + cf^ are degenerate with each other at the one- 
loop order since we only track the growing modes, and degenerate with other parameters at 
higher loop order. While the shear viscosity parameter cj.^ affects the vorticity, which is a 
somewhat small effect. As analyzed in Ref. |32] one can directly evaluate the stress-tensor 
from the microphysical theory, i.e., from N-body simulations to extract such parameters. For 
smoothing scale A = 1/3 [h/Mpc] at z = it is found + fc^ ^ 9 x 10~^c^, (/ is the logarithmic 
derivative of the growth function). This direct measurement can also be obtained from matching 
to the power spectrum at some renormalization scale, resulting in a consistent value and a positive 
check on the validity of the theory. 

We establish recursion relations for the density fluctuations and velocity field, allowing us 
to insert this measured value of the fluid parameter and obtain correlation functions. These 
parameters carry A dependence which balances the A dependence of the cutoff on the loops. If 
we take A to large values the fluid parameters approach a finite quantity, representing the finite 
error made in the standard perturbation theory. In particular the fluid parameters provide the 
following relative correction to the power spectrum (suppressing the time dependence and C(l) 
factors) 

where PL{k) is the linear power spectrum. Note that the pressure and viscosity act together 
to reduce the power spectrum by acting oppositely to gravity. This simple, but entirely real 
and rigorous correction to the power spectrum is essential to explain the observed shape of the 
baryon-acoustic-oscillations in the power spectrum relative to standard theory. 

The outline of the paper is the following: In Section [2] we describe the basic theoretical 
setup. Operating in the Newtonian approximation in an expanding universe, we smooth the 
Boltzmann equation to obtain an effective fluid for cold dark matter. We recapitulate how to 
match its parameters to the microphysical results from N-body simulations. In Section [3] we 
solve the theory recursively for a matter dominated universe by extracting the time dependence. 
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and lift this to ACDM in an approximate way also. This allows us to semi-analytically derive 
the power spectrum at one-loop order. In Section |4] we present our numerical results for the 
power spectrum and compare to linear theory and standard perturbation theory. In Section [5] 
we discuss the fluid's parameters, corrections from collisions, wave-like behavior, higher order 
moments, and the velocity field. Finally, in Section E] we summarize the effective field theory 
and its role in cosmology. 

2 Effective Fluid 

2.1 Newtonian Approximation 

Cosmological perturbation theory around a fiat FRW background may be performed in many 
gauges. One example is the Newtonian gauge. Scalar modes are captured by the following form 
for the metric 

ds^ = -dt^{c^ - 20(x, t)) + a{tf{l - 27/;(x, t)/c^)dx^ (3) 

were a is the scale factor and x is a co-moving co-ordinate. The Newtonian approximation is a 
valid description for non-relativistic matter in an expanding background on sub-horizon scales, 
and will be sufficient for our purposes as we will study evolution of matter after equality. In this 
limit, only plays a role and not ip. Here is the Newtonian potential, sourced by fluctuations 
in matter density [33] 

V20 = 47rGa2(p(x,t)-pfe(t)) (4) 

where p is the matter density, which is a combination of dark matter and baryonic matter, and 
Pb = (p) is the background value, with pb{t) oc l/a{t)^. The Hubble parameter is determined by 
the Friedmann equation 

H{tf = ^iptit)+p,ac) (5) 

in a flat ACDM universe. Here we allow for vacuum energy in pvac, which we assume to be the 
cosmological constant of general relativity, as is consistent with all current data. The ACDM 
concordance model indicates that this is a valid description of the universe for all times well 
after matter-radiation equality. Furthermore, after the baryons decouple from the photons, the 
baryons tend to trace the dark matter on large scales, leading to a single non-relativistic fluid 
that we will describe. 

2.2 Phase Space Evolution 

We treat dark matter as classical point particles and ignore its quantum nature. This is a very 
good approximation for most dark matter candidates, but can breakdown for extremely light 
axions which organize into a state of very high occupation number, with quantum pressure and 
associated sound speed Cg ~ [31]. For QCD axions, this is ignorable on large scales (since 
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it vanishes for small k) and will be ignored here; see Section 15.41 for further discussion. Indeed 
N-body simulations are ordinarily done with a set of classical point particles. At each moment 
in time the output is a set of N- vectors which we label n, with co- moving co-ordinates x„, and 
proper peculiar velocity v„. 

Let /„(x, p) be the single particle phase space density defined such that /„(x, p) (i'^xc/^p is 
the probability of particle n occupying an infinitesimal phase space element. For a point particle, 
the phase space density is 

/„(x,p) =(5|)(x-x„)(5|)(p-mav„) (6) 

(where both x and p are co- moving). By summing over n, we define the total phase space 
density /, mass density p, momentum density vr*, kinetic-tensor cr*-' as 

/(x,p) = ^5|)(x-x„)(5|5(p-mav„) (7) 

n 

p(x) = ma-^j rf^p/(x,p) 

= ^ma-^5|,(x-x„) (8) 

n 

7r^(x) = a~' [ rf3pp7(x,p) 



5^ ma 3<5|,(x-x„) (9) 



(j*^'(x) = m '^a ^ J (i^pp'>V(x, p) 

= ^ma-3<<4(x-x„) (10) 

n 

The Newtonian potential is sensitive to an infrared quadratic divergence in an infinite homo- 
geneous universe. To isolate this divergence we introduce an exponential infrared regulator with 
cutoff n (a "mass" term) and will take the /i — )■ limit whenever it is allowed. The Newtonian 
potential (p is 



J |x - x'l 



ax — x„ 



2 / ^3„/ P(xO - Pfe ^-^lx-x'l 



(x) = -Ga' d'^ "^r e 
J |x-x'| 

= l^^n^ -T— (12) 



Note that 



fx^ ^^0n — ^ —4:71 Go? ph as /i — )■ (13) 



The /c-space version of the Newtonian potential is 

a(A;'' + /i^) 

<^(k) = ^0„(k) + l^^(2vr)35|,(k) (15) 

jjj 

n 

where the final term evidently subtracts out the zero-mode. 
2.2.1 Boltzmann equation 

Cold dark matter candidates, such as WIMPs and axions, have very small scattering cross 
sections with itself and standard model particles. Here we make the approximation that we can 
ignore the scattering altogether (see Section [5^ for further discussion). Restricting our attention 
then to coUisionless classical particles interacting only via gravity, the Boltzmann equation is 

p"^ + r>v|:)/,. = o (16) 

In the Newtonian limit in a fiat FRW expanding background, the coUisionless Boltzmann equa- 
tion becomes 

D/^ ^ 9/^ _P_ ^ Q^N 

Dt dt ma? ' 9x I ^ (9x j ' (9p ^ ' 

where the self-force has been subtracted out of the sum over n in the last term. By summing 
over n we have 

o = t;7 = 7^ + — 2 ' IT ' ir~ 1^ 

Dt at ma'^ ax ^ ax ap 
where the final term now involves a double summation over n and n. 



2.3 Smoothing 

From this we would like to construct various quantities that describe an effective fluid at large 
length scales [31]. Most importantly, we need the effective stress-tensor r*-' that is sourced by 
the short wavelength modes. In order to deflne a fluid we must perform a smoothing of the 
output data from the simulation. 

To this end let us deflne the following Gaussian smoothing function 



8 



(here A is a smoothing scale and should not be confused with the cosmological constant) which 
is normalized such that J (i'^xiy(x) = 1. Of course our results will not depend on the choice of 
smoothing function, but the Gaussian is chosen for convenience. (In fact we will later include 
numerical results for the sine function corresponding to a step function in fc-space). In /c-space 
the Gaussian smoothing function is 



Wj,{k) = exp(^-^^ (20) 

This will smooth the fluid quantities on (co- moving) length scales ^ A"^, acting as a cutoff 
on modes k > A. We should choose A < kj^^, were /c^l is the wavenumber where modes have 
become non-linear, so that we integrate out the non-linear modes. Given the form of W\, we 
can estimate a rough value for the smoothing scale in position space as Auv ~ 27r/A. 
For certain observables 0(x), we will define the smoothed value by the convolution 



OKx) = [0]a(x) = j dVW^A(x-x')0(x') 



(21) 



The smoothed versions of the phase space density fi, mass density pi, momentum density tt^, 
stress-tensor af, derivative of Newtonian potential di(pi are 

/Kx,p) = ^iyA(x-x„)5|5(p-mav„) (22) 

n 

pKx) = ^ma-3l^A(x-x„) (23) 

n 

TT^x) = ^ m < Wa{^ - xj (24) 

n 

(x) = ^ m a-=^ < vi W^i^ - x„) (25) 

n 

The "1" substrict indicates that these only depend on the long modes. Similarly, the smoothed 
version of the Newtonian potential (pi is 



AlX — X»; 



V2 



e-M|x-x„| ^26) 



mG 

0«,n(x) = ■ rErf 

a|x - x„| 

n 

Here we used the Gaussian smoothing function to explicitly evaluate 0/(x) in terms of the error 
function Erf (2;) = ^ dt e~* . 
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We now write down the smoothed version of (fTSl) by multiplying it by Wa and integrating 
over space 







Dt 



dfi ^ P dfi 
dt m 9x 



(2^ 



where we used linearity to express the first two terms directly in terms of //. However the third 
term is more complicated. The equations of motion are obtained by taking moments of this 
smoothed Boltzmann equation, i.e., 







21 

Dt 



(29) 



J A 



One often subtracts out traces of the tensor structure of the higher order moments for con- 
venience. Here we will only make use of the zeroth and first moments, so this detail is not 
important for us. 



2.4 Effective Continuity Equation 

The zeroth moment gives the continuity equation, which we find to be 

Pi + 'iHpi + -di{piv\) = Q 
a 



Here we introduced the velocity field 



PKx) EnW^A(x-X„) 



(30) 



(31) 



The continuity equation relates the zeroth moment of the phase space distribution, p^, to the 
first moment of the phase space distribution, piv\. We now turn to the next moment of the 
Boltzmann equation to obtain an equation for the velocity field itself. 

2.5 Effective Euler Equation 

The first moment gives the Euler equation, which we find to be 



a ' a 



a Pi 



where the effective stress-tensor that is sourced by the short modes is given by 



(32) 



(33) 
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The Euler equation relates the first moment of the phase space distribution vl, to the second 
moment of the phase space distribution, a^-' . Here is a type of kinetic dispersion and is 
a type of gravitational dispersion, namely 

«f = - pivivi (34) 
" SvrGa^ + SyrGa^ ^^^^ 

where 

= y rfVWA(x-x) [9i/0(x ) 9jv0(x ) 

-5^9,,0„(x')9,v0„(x')] (36) 

n 

Note that we have subtracted out the self term in w\\ and used = 4:TiGa?{p — pb) and 
V^0/ = 4:TTGa?{pi — Pb) to express $/ in terms of and (pi. In the limit in which there are no 
short modes, it is simple to see from the definition of and that they vanish in this limit. 
More on the mathematical details of this are given in Appendix |Al 



2.5.1 Derivative of Stress- Tensor 

In the general relativistic theory the absolute value of the stress-tensor will act as a source for 
gravity, and could, in principle, be important at the scale of the horizon. This includes the trace 
of the stress-tensor that is a type of pressure, which we discuss in Appendix [B] On sub-horizon 
scales, however, a non-relativistic analysis is applicable in which the right hand side of the Euler 
equation only involves the derivative of the stress-tensor: 

j; = -i-9,nA (37) 
a Pb 

here we divided by the background density pb for convenience. This quantity will be quite 
important in out analysis, and takes on the following explicit form 

a Pb Jl = dj{al^ - piv\vl) + ^ [p„ 9i0fJ^ - pidi(t)i (38) 

The smoothed quadratic form can be expressed as 

[Pn d4n]t, = ~ "i^ ( 1 + /^|xn - x^|)e-^l— ''lTyA(x - x„) (39) 

^ I x^, x„ I 

which requires one to perform a double summation over n and n, which can be computationally 
expensive. Notice, however, that while the stress-tensor involves an integral over the gravita- 
tional potential in eq. f p6|) . the derivative of the stress-tensor does not require this (by making 
use of the Poisson equation and then integrating over the delta-functions). 
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2.6 Fluid Parameters 



The effective stress-tensor [t^-']k comes from smoothing over the short modes and therefore is 
"sourced" by the short modes. However, the important quantities that will arise when we later 
compute A^-point functions, such as the 2-point function, involves correlation functions of the 
stress-tensor with the long modes. 



and 5i (x) 



pb 



1 



(40) 



(note that 0; is determined as a constrained variable through the Poisson equation V^0i = 
47rGa^PbSi). For instance, the expectation value {[t^^a) is some background pressure. More 
importantly though is the mode-mode coupling. Coupling between long modes is connected to 
the non-linear terms in the continuity and Euler equations. While coupling between long and 
short modes is connected to the stress-tensor, which generates non-zero correlation functions 
{[t^-^]a^i) and ([t^-'Ja "wf ). The long-long and short-short mode couplings are represented by 
vertices in Fig. [H Note that the dark matter gas is non-thermal, and this indicates that the 
stress-tensor cannot be derived by some analytical thermal argument, such as would be the case 
for air under ordinary conditions. However, the reason one can make progress is to recognize 
that the stress-tensor organizes itself in terms of length scales with coefficients that come from 
matching; this organization is guaranteed by the principles of effective field theory. 

In order to extract this dependence of the stress-tensor on the long modes, we implicitly write 
the stress-tensor as an expansion in terms of the long fields, whose coefficients are determined 
by various correlation functions. This will involves a type of pressure perturbation term oc S^^ Si, 
a shear viscosity term oc djvl + divj — dkVi, and a bulk viscosity term oc S^^ d^vf. Demanding 
rotational symmetry, we write a type of effective field theory expansion for the stress-tensor as 



pb 



Ha 



(41) 



where 7 would correspond to the ratio of specific heats in an ordinary fiuid (e.g., 7 = 5/3 for 
an ideal monatonic gas) but here it just parameterizes the background pressure term, Cg is a 
sound speed, and Cgv , c^v are viscosity coefficients with units of speed. Note that Cg, Csv, Cbv are 
allowed to depend on time, but not space. Our fiuid coefficients are related to the conventional 
fluid quantities: background pressure pb, pressure perturbation 6p, shear viscosity 77, and bulk 
viscosity ( by 



Pb 



clpb 



1 

AH 



6p = clpb^u 
^ ~ H 



(42) 
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Figure 1: Vertex for interaction between long-long mode coupling or long-short mode coupling. 

In addition to what is included in (HTl) there is an entire tower of higher order corrections 
carrying the appropriate rotational symmetry, guaranteed to exist by the principles of effective 
field theory. These will be parametrically suppressed at low wave number k compared to the 
non-linear wavenumber k^i, and will not enter to the order we shall work (which will be 
and so shall be ignored in the present discussion. Here Ar*-^ represents stochastic fluctuations 
due to fluctuations in the short modes, with (Ar*-') = 0. We will return to this issue in Section 

El 

For convenience, let us define the following quantities from the stress-tensor 

j; - (43) 

4- s -i-9,9,.[r-']_, = 4J,7a (44) 

0, Pi) 

Ai = -^d,d,[T'^]^ = dai/a (45) 
Now we introduce the dimensionless velocity divergence 
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Then according to our ansatz (14T|) (and ignoring stochastic fluctuations for now) we have 

a Ji = cl dA + \cl a,ef + + cL) d,Qi (48) 



a^Af = cl dM + -cl dud^Qf +[^ + cl] dAQi (49) 



T 

a^Ai = cl dHi + (4 + cl,)d^Qi (50) 
a^Bi = cld'Qi, (51) 

In order to extract the coefficients, we multiply each of these by the functions on the right 
hand side and then form a position space correlation function (...), say (V'i(x + x')?/^2(x')). We 
will need the following set of correlation functions 

(52) 
(53) 
(54) 
(55) 
(56) 
(57) 
(58) 
(59) 

By rearranging, we find the following expressions for the fluid parameters 



Pa5\ 


[X) 


= (A,(x4 


-xO^KxO) 


PAe\ 


[X) 


= (AKx4 


-x')eKx')) 




[X) 


= (4^(x 


+ x')ef(x')> 




[X) 


= (fi^(x^ 


-xOe^x')) 


P55\ 




= (5/(x + 


xO^Kx)) 


Pse^ 




= (5Kx + 


x')0Kx')) 


Peel 




= (eKx^ 


-x')eKx')) 






- (ef (X 


+ x')er(x')> 



C2 



2 ^ PAejx) d^Psejx) - Pas{x) d^Peejx) 

{d^Pse{x))ya^-d^Pss{x)d^Pee{x)/a^ ^ ' 

2 ^ Pas{x) d'Psejx) - PAejx) d'Pssjx) 

{d^Pse{xW/a^ -d^Pss{x)d^Pee{x)/a'' ^ ^ 

4 Py^k^Qk^{x) - PAejx) 
3 ^^Pek^ek^{x)/a^ - d^Pee{x)/a^ 

Ppeix) 
d^Peeix)/a^ 



(62) 



where 



is the sum of the viscosity coefficients. The final result for each of the fluid coefficients will surely 
have some spatial dependence, so one should take the x ^ A^l limit of the final result to extract 
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a constant value. Note that in eq. (!62|) we have provided two alternate expressions for c^„. In 
the linear theory, 5 ~ G, allowing one to approximate the sum 



2 2 _ PAejx) _ PAsjx) 

''^ 52Pee(x)/a2 ^ d^Ps6{x)/a?- ^ ^ 

These expressions may also be given in /c-space by taking the k <^ k^^ limit. To do so we 
Fourier transform each of the correlation functions. We define the Fourier transform as 

0(k) = y d3xe-'''-"0(x) (65) 

By translational invariance, each /c-space correlation functions takes the form 

(V'i(k)V2(k')) = (27r)35|,(k + k')P^,^,(A;) (66) 

where P^^ (k) is the Fourier transform of P^j (x). This allows us to write the fluid parameters 
as 

2 ^ PAsik) Psejk) - PAsjk) Peejk) 

-k^Pse{kY/a^ + k^Pss{k)Pee{k)/a? ^ ^ 

2 ^ PAs{k)Pse{k)-PAe{k)Pss{k) 



C2 



-PP,e(A:)V«' + k'Pssik) Pee{k)/a^ 
4 PAfc.0fc»(A;) - Pyie(fe) 

3 -k'^PeMQ^.{k)/a? + k'^Pee{k)/a? 

Pseik) 
-PPee(A;)/a2 



(69) 



We shall soon who that the combination + cl is most important for the leading correction to 
the power spectrum. Whether it is possible to measure each of these parameters individually, or 
only certain combinations, is discussed later in Section 15.21 



2.7 Relative Size of Terms 

Let us recall the primordial power spectra. In the first part of this subsection we will focus on 
a pure Einstein de Sitter universe, i.e., ignore the turn-over in the power spectrum due to the 
transfer function T{k) (see next Section for its description). We will then include comments on 
the change that occurs when the transfer function is included, and all our numerical results in 
the latter part of this paper will be for the real universe including the full transfer function. 

The primordial power spectrum in the Newtonian potential is approximately scale invariant 
(in this section we will suppress factors of 27i) 

P^ik) - -^r- (70) 
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For sub-horizon modes that entered the horizon in a matter dominated era, the Poisson 



equation gives = —\H'^a? 5k and the corresponding power spectra is 



Pso(k) - pi) 
This means that the characteristic value of 5i in position space, on a scale set by /c, is 

5i ~ ^/WlUk) ~ 10-^ (^-^^ ' (72) 

This estimate is valid for k < keq- For k > k^q the rise in 5i is only logarithmic. 

Let us now compare the relative size of the terms that appear in the Euler equation. We use 
the Hubble friction term Hv] as the quantity to compare to. The relative size of the pressure 
c^diSi/a or viscosity c^did^Vi / {Ha?) ~ c^di5i/a is 



Pressure, Viscosity ^, k5i/a 

rsj ■ 

Hubble Friction Hvi 



(73) 



k V cl 



_ IJ^\ ^ ^2:^ 5 (74) 



The relative size of the non-linear piece of the velocity convective derivative v\djv\/a is 

non-linear Velocity kvf/a 



Hubble Friction Hvi 



(75) 



ck \ ^ 



~10-(-j~*, (76) 

Note that since we expect ~ ~ 10~^c^ for a pure Einstein de Sitter universe, then the 
pressure, viscosity, and non-linear velocity piece appear to be comparable on all /c-scales and 
all are ~ Si. This also means that the terms 5^^\J^^^^ and J^^'^\ which appear in the Taylor 
expanded Euler equation, are suppressed by a factor of ~ 6f and so are all higher order again. 
These estimates are consistent with the fact that we should only probe scales larger than the 
non-linear scale, i.e., k~^ > k~^ ~ 10~^'^^ c/ (Ha). 

On the other hand, our analysis is not applicable in the regime approaching the horizon scale 
where general relativistic (GR) corrections will be important. For instance, GR will change 
the partial derivatives to covariant derivatives leading to Hubble corrections. This leads us to 
estimate 

GR Correction f HaV 10"^ 



- ~ -— (77) 
Newtonian Approximation \ck J 6i 

In other words, we should only probe scales smaller than Hubble, i.e., k~^ < kjj^ ~ c/{Ha). 
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For baryon-acoustic-oscillations the relevant scale is roughly k^^^ ~ -^{Hdecddec) ^ at the 
time of decoupling, where is the sound speed of the photon-baryon plasma (the associated 
peak in the fluctuations will be at a somewhat smaller scale when all modes are properly in- 
cluded). Red-shifting to today gives an estimate for 5i that is somewhat less than 1 on this scale. 
So this appears to fit nicely in the window where our approximations are valid. However, this 
estimate is overly simplistic, as it ignores the turn-over in the power spectrum for modes that en- 
ter during the radiation dominated era as described by the transfer function T{k). Nevertheless, 
this qualitatively sets the basic hierarchy. 

Turning then to the various length scales in the real universe: the characteristic length scale 
for baryon-acoustic-oscillations is Abao ~ 120 Mpc (see Fig. Hj). The non-linear scale is not 
sharply defined, since the power spectrum turns over to a logarithm for modes that entered in 
the radiation dominated era, but a characteristic value is A^l ~ 10 Mpc. This suggests we put 
the cutoff scale on the effective fluid at Auv ~ 2ti / K ~ 20 Mpc, or so, in order for it to be just 
above the non-linear scale. The Hubble scale du ~ 4Gpc is a lot larger and so too is the scale 
of equality Agg ~ 600 Mpc. A summary of this hierarchy is 

Anl < Auv < Abao < ^eq < du (78) 

Although one should be careful here; although the baryon-acoustic-oscillation scale is ~ 120 Mpc, 
its width is much smaller ~ 20 Mpc. So one may actually need a somewhat smaller Auv than 
~ 20 Mpc to fully resolve the baryon-acoustic-oscillation peak. But this then starts to push up 
against the non-linear scale ~ 10 Mpc. However, we will eventually send Ayy — )■ 0, once we have 
removed the cutoff dependence, so this is not necessarily a problem. 

3 Perturbation Theory 
3.1 Linear Power Spectrum 

Inflation generates the primordial power spectrum, which we assume to be Gaussian 

(5i(k)5i(k')) = (27r)35|,(k + \^)PUk) (79) 

where the "L" subscript indicates that the initial fluctuations are in the linear regime. Here 
nf oc fc"^, with Us — 1 for a scale invariant spectrum. After inflation, one draws 5(k) from 
Pinf(A;), as well as photons, electron fields etc, and evolves with a program such as CMBFAST 
or CAMB to sometime after recombination. This should be adequately captured by the linear 
evolution, but should be a fully relativistic calculation. 

Lets call the initial scale factor after inflation Oj and the late time to which we evolve under 
linear evolution by aiate- The final density fluctuation will be related to the initial density 
fluctuation by the transfer function T{k) and the growth function D{a). One defines the transfer 
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function as the ratio of the gravitational potentials on a given scale to that on the large scales, 
i.e., 

<Pls[k,aiate) 

In particular for modes that enter during the matter dominated era, we have T{k < keg) ~ 1. 
The transfer function decreases for modes that entered in the radiation dominated era. In 
particular 



with wavenumber at equality [35] 



keg ^ 0.073 Mpc"^ oh^ (82) 



Plus, there are corrections from baryon-acoustic-oscillations etc on the spectrum which CMB- 
FAST or CAMB should provide. 

At the linear level, we must then simply multiply by the growth function D{a) to reach 
today's spectrum. We express the evolution in the Newtonian potential as follows 

Mk,a) = ^<PUk)nk)^ (83) 

where the 9/10 prefactor is from the Sachs- Wolfe effect, (p^^i is the primordial fluctuation, T{k) 
is the transfer function, and D{a) is the growth factor normalized to the scale factor a for 
convenience. Later we will generalize the growth factor from D{a) — > D{k, a) to account for the 
/c-dependence in the resummed linear theory, but lets suppress that for now. 

For sub-horizon modes, the Poisson equation gives the following relationship between density 
and Newtonian potential 

_ 2P0La , . 

- iTFf 772 1^4) 

where flrn,o and Hq are today's values. So the density fluctuation is given in terms of the 
primordial fluctuations, transfer function, and growth factor by 

^Lik, a) = ^^f^<p,^,{k)T{k)D{a) (85) 

The primordial power spectrum generated during inflation is 

For a power law, we write the power spectrum as 




lOL 

10" 



0.001 0.01 0.1 

Wavenumber k [/i/Mpc] 



Figure 2: Linear power spectrum of density fluctuations Piik) computed from CAMB, with 
rig = 0.96, 2 = 0, Qm = 0.226, flk = 0. The plot for k < keg shows the approximate scale 
invariance of the spectrum. 



where we followed [35] in the definition of the amplitude, denoted Sh- Combining the above, we 
have the linear power spectrum in the density as 



PL{k,a) 



0^2 r2 



7Tns+3 



T\k) 



Dia] 



Da 



From running CAMB, a plot of this linear power spectrum is given in Fig. |2l The dimensionless 
variance is defined as 

A^fc, a) ^ ^'^^^^ (89) 



Note that on Hubbies scales today; A^(ifo,fl 



27r2 

1) = 6fj, which explains the unconventional 



normalization chosen in P^{k). The measured value of the amplitude of density fluctuations on 
the scale of the horizon is ~ 1.9 x 10~^ [36] . 

The standard deviation ^^{k^a) is a measure of the fluctuations in 5 on a scale fc, which 
we plot in Fig. [3J Since these fluctuations become larger than 1 at high k, the theory is non- 
linear in this regime, which sets a non-linear scale of 

/cnl ~ 0.5 [h/Mpc], or so. This defines 
the perturbative regime in which the effective fluid description is applicable for k < /c^l and the 
non-perturbative regime in which the effective fluid description breaks down for k > fc^L. The 
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Wavenumber k [/j/Mpc] 

Figure 3: Linear standard deviation of density fluctuations ^5{k) computed from CAMB, with 
Hs = 0.96, z = 0, Qm = 0.226, Qk = 0. The plot indicates that the evolution is perturbative for 
small k and non-perturbative for high k. 



scale /cnl acts a type of "coupling" in the effective field theory and organizes the expansion into 
powers of k/k^^. 

By reverting back to position space we can define the correlation function ^. Using statistical 
homogeneity and isotropy it is related to the variance by 

e(r,a) = (5(x,a)5(x + r,a)) = /"rfln A; ^^^^A?(A;, a) (90) 

J kr 

and is in Fig. HI The correlation function evidently includes the baryon-acoustic-oscillation peak. 
Its precise shape is subject to non-linearities that we would like to compute accurately; although 
this is analyzed most cleanly in fc-space, which will be our focus. 



3.2 Evolution Equations 

In order to study non-linear corrections, we begin by recalling here our equation of motion for 
the velocity field for the stress-tensor ansatz we made earlier 

v\ + Hv\ + vid,v\ + -Ui = -\cl OA + ^,d'vi + ^^^d,d,vi - A (91) 
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Distance r [Mpc/h] 

Figure 4: Linear density correlation function ^(r) computed from CAMB in a ACDM universe, 
with Us = 0.96, z = 0, Qm = 0.226, Qk = 0. This clearly shows the baryon-acoustic-oscillation 
peak at r ~ 120 [Mpc/h]. 

where A J* = p'^^djAr'^^ /a. Which is complimented by the Poisson and continuity equations. 
Our ansatz for [t''^]a suggests that the right hand side of (PT]) should be multiplied by an overall 
pref actor 1/(1 + 5^). But since we have only included linear terms in our ansatz for [t'^-']a it would 
not be self-consistent to include this prefactor. Furthermore, the leading order from such terms 
would enter parametrically at 3rd order in an expansion in powers of the linear density field 
5^^^ as c^V^fJ^^-* 5^^\ This would correct the two-point correlation function as (c^V^^^^^ (5*^^^ 5^^^), 
which obviously vanishes when the primordial spectrum is Gaussian (or just even in (5(1)). Hence 
we can drop such corrections. 

3.2.1 Curl of Velocity 

Before examining the density fluctuations in detail, let us briefly mention the vorticity. The curl, 
or vorticity, of the velocity field 

w/ = V X Vi or in /c-space = i k x v/ (92) 

is determined by taking the curl of the Euler equation. We use the vector identity 

V X (v, ■ V)vO = -V X (v; X (V X V/)) (93) 
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to obtain the non-linear vorticity equation 



Let's first discuss this at the hnear level w*^^) where we ignore the right hand side. Even in 
the absence of viscosity, the vorticity w^^^ is being driven to zero in an expanding universe as 
w/ oc 1/a, which is a well known result. In the presence of viscosity, this happens all the more 
rapidly (assuming c^^ > 0). The curl of velocity in a matter dominated universe is plotted in 
Figure |3 This means that at the linear level, studied at late times, we cannot see the effect of 
the shear-viscosity as the transient vorticity will have decayed away. The non-linear term on the 
right hand side also vanishes when = and therefore vorticity is not generated, unless it is 
present initially (although the curl of the stochastic fluctuations V x A J could alter this). 

However, other non-linear terms that we have neglected, such as the overall prefactor 1/(1 -|- 
6i), will generate vorticity. These are required for the shear viscosity to have an explicit effect, if 
we cannot track the initial transient viscosity. On the other hand, the combination of viscosities 
cl = c^^ + cl^, will appear in the divergence of the Euler equation (see eq. f l98|) ) and can play an 
important role at late times. For the present discussion, let us assume the vorticity is negligible. 
This is consistent with simulations and observations which find the vorticity to be small, albeit 
non-zero. 



3.2.2 Divergence of Velocity 

The divergence of the velocity 

6i = V ■ or in fc-space 6i = ik ■ vi (95) 

is coupled to the density fluctuation 6i. Lets write down the coupled equations using conformal 
time and the associated Hubble parameter 

n=^-^ (96) 



a{t) a dr 

In the absence of vorticity, and ignoring stochastic fluctuations for now (see Section [D] for its 
inclusion) the evolution equations for the pair 6i, 9i are found to be 

^ + ^^ = -J ^3^(^^ k')5,(k - k')^,(kO (97) 

^ + Wi + ^H'Qm^i = - / 7^/3(k, k')ei{k - k%ik') + clkHi - ^-^01 (98) 
dr 2 J [Zir)'^ ti 
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Figure 5: Curl of velocity w (normalized to some initial value) in a /c-mode as a function of time 
in a matter dominated era in the linear approximation. Blue is for c^^ > with re-summation 
and red is for cf„ = 0. This shows the decay in vorticity over time at this order of analysis. 



which comes from taking a divergence of the Euler equation. Here 

kk' 



a(k,k') 
/3(k,k') 



e k'(k 



2|k'P|k-k 



/|2 



(99) 
(100) 



appear as the kernels of the above convolution integrals. The convolution integrals on the right 
hand side of eqs. ([HZHM]) arise from long-long mode coupling, while the Cg and Cy terms arise 
from long-short mode coupling. 



3.3 Recursion Relations 

Since we will be evolving the universe fully linearly during the radiation dominated era, until 
after decoupling, the subsequent evolution will be during a matter dominated era. Although the 
late time behavior with a cosmological constant will alter this quantitatively, let us ignore this 
for now and study the matter dominated era for analytical simplicity; the generalization to the 
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ACDM case will be performed in Section 13.61 In the matter dominated case we have 



2 (r^^ 



fi„ = 1, n = -, a= (^-J (101) 

Given this, we would like to extract the time dependence in the problem, by performing a self 
consistent expansion in the scale factor, as Hubble is the only time scale in the problem. 

Recently in Ref. [32] we exploited the use of Green's functions to capture the time dependence, 
which is a very powerful technique. This led to result in terms of integrals over factors of Green's 
functions, etc. Here we would like to extract the time dependence in a more explicit and intuitive 
way to compliment the previous powerful results. Since the density field and velocity field is 
small at early times, we expand our fields in powers of the scale factor as follows 

oo 

5Kk,r) = J]a"(r)5„(k) (102) 

n=l 

oo 

^,(k,r) = -?/^a"(r)^„(k) (103) 

where we have suppressed "1" subscripts in the perturbed expansion on the RHS, although they 
are all long-modes. This cleanly separates out the time dependence and the /c-dependence. 
Ordinarily, the higher order contributions to the density fields (5„ can be counted in powers of 
the linearized density fluctuation 5i, namely (5„ ~ ~ 5". However, the presence of the fluid 
corrections alters this simple counting here. Note that the leading order terms satisfy 5i = 6i, 
with growth factor D{a) = a. The linear density fluctuation 5i(k) is therefore drawn from the 
time independent power spectrum 



Pnik) ^ ^^^W^{kf (104) 



J^ns 



= 2n'6l--^,T\k)W4kr (105) 

where the factor WA^k)"^ comes from smoothing and T{k) is the transfer function that includes 
the full effects from the radiation dominated era as computed by a program such as CMBFAST 
or GAME. 

If we substitute the expansions for 6i and 6i into the continuity and Euler equations, we can 
equate powers of the scale factor a, giving a pair of recursion relations for 6n and 9n- This will 
depend on the time dependence of the fluid parameters Cg and c^. In a matter dominated era, 
we will show in the Section 13.51 that the approximate time dependence of the fluid parameters 
is that they increase with the scale factor a. To capture this we introduce 



c 



2^a, A) = a<o(A), cl{a,A) = a<o(A) (106) 
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and we define the following time independent dimensionless parameters 



„2 h.2 2 1.2 



-1/2 ' 
rtn 



(107) 



where the A dependence is implied. For n > 1 we find the following set of relationships between 
the fields at different orders 

yl„(k) = n5„(k)-^„(k) (108) 
5„(k) = 35„(k) - {2n + l)9M - 2C,,o(A:)5„-2(k) - 2C„,o(A;)^„_2(k) (109) 

where the left hand side are the following non-linear integrals 



A„(k) = J^sJ ^'^2 5|>(ki + k2 - k) 

n-l 

a(k, ki) ^ 9miki)6n-mi^2) 



;iio) 



m=l 



(2^ 



d=^A;2 5|5(ki + k2-k) 



2/3(k,ki)5];^^(ki)^,,_„(k2) 



:iiii 



m=l 



These relationships allow us to express the m = n^^ value of the fields in terms of the m < n^^ 
value of the fields. Namely, we have the following recursion relations for > 1 



5n(k) 



1 



, (2n+ l)AJk) - BJk) 
(2n + 3)(n- 1) ; nv ; nv ; 

-2Cs,o{k)6n-2{k) - 2a,,o(A;)^„_2(k) 

^- f3A„(k) - nBJk) 

(2n + 3)(n- 1) L ^ ^ 

-2nCs,oik)5n-2{k) - 2na,o(A;)^n-2(k) 



(112) 



:il3) 



with starting values ^i(k) = 6i{k). This allows us to in principle solve for all the higher order 
fields in terms of 5i(k). 

The solution for Sn and On. can be expressed in terms of kernels Fnj and Gnj, rather than 
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left in terms of the stochastic variable 6i. Lets write our fields as 



5„(k) = X:y"^...y"rf'g,4(qi + --- + q.-k) 



(2vr)3--- 
F„j(qi,...,qj)5i(qi)...5i(qj-; 



^n(k) = |j|A^...|rf3^,.5|,(q, + ... + q,.-k) 
G„j(qi,...,qj)(5i(qi)...(5i(qj) 



(114) 



;ii5) 



Then by substituting into flll2H113"l) we establish recursion relations for F„j and Gn.j with 
starting values Fi^i = Gi^i = 1. At the one- loop order, we will find that only Gn,n and 
Gn,\ will enter, if the primordial fluctuations are Gaussian. So we report on their values here. 
Firstly, and G„,.„ are independent of the fluid parameters and are given recursively by 



iqi, 



n-1 

^ (2n + 3)(n-l) 

m=l ^ ' ^ ' 

[(2n + l)a(k, \^\)Fr,-m,n-mkPim+\, 

+2/3(ki,k2)G 

n— m,?i— ml,Qr7i+l ) • • • ; Q.n 



:ii6) 



n-1 



Gn,n{(il, ■ ■ ■ , Qn) — ^ 



qr: 



(2n + 3)(n-l) 
[3a(k,k 

l) -^n— m.n— m(qr?i+l ) • • • ; qn) 

+2n/3(ki,k2)G 



:ii7) 



On the other hand, Fn^i{k) and Gn,i{k) are determined entirely by the fluid parameters. We 
find them to be the following products {n > 1) 



F,,j(k) 



n 

m=3.5,... 



-2(C,(A;) + (m-2)a(fc)) 
(2m + 3)(m-l) 

-2(C,(fc) + (m-2)a,o(fc)) 
(2m + 3)(m-l) 

and we find = Gn,i = for n even. Note that Gn,i = n Fn^i. 



m=3,5,. 



for n odd 



for n odd 



'118) 



:ii9) 
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3.4 Power Spectrum 

We will go to one-loop order in the power spectrum. This will require the second and third order 
density corrections. At second order, we have the following (symmetrized) kernels 



2,2 ' 



ki,k2 



G2,i{k) 



5 ^ 2 (ki-k2)^ ^ ki-k2 

3 ^ 4 (ki ■ k2)^ ^ ki ■ k2 

7 7 k-^ k<2 2 






~^ ^2 

"^1 "-2 
1 1 

Av]^ ft.2 



(120) 

(121) 

(122) 
(123) 



At third order we need F3 3 , 6*3^3, -F3 1, and G^^i which we find to be the following (unsym- 
metrized) kernels 



^3,3(qi,q2,q3) 



7a(k, qi)F2,2(q2, qs) + 2/3(qi, q2 + q3)G'2,2(q2, qs) 



+ (7a(k, qi + q2) + 2/3(qi + q2, q3)G2,2(qi, q2) 



G'3,3(qi,q2,q3) = — 3a(k,qi)F2,2(q2,q3) + 6/3(qi,q2 + q3)G'2,2(q2,q3) 
io L 



-(3a(k, qi + q2) + 6/3(qi + q2, q3)G2,2(qi, q2) 
1 



(C,,o(A;) + a,o(A;)) 



The two-point function for 61 defines the smoothed power as follows 



(124) 

(125) 
(126) 
(127) 

(128) 



We now substitute in the expansion fll02p for 61 in powers of the scale factor to give the two-point 
function the form 



{6i{k,T)Si{k',T)) = a^(r)(5i(k)5i(k)) + 2a3(r)(5i(k)52(k')) 
+ a^r) [2(5i(k)53(k')) + (52(k)52(k'))] + • 



(129) 



In order to organize this into an expansion in power spectra, let us define the power spectrum 
that contributes at n'^^ order as 



(5„(k)5„_^(k')) = (27r)35^(k + k')P^.-„(A;,r) 



(130) 
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By inserting this into f ll29p we obtain the following expansion for the power spectrum 



+ a\T)[2Pr,{k) + P22{k)] + ... (131) 

We assume that the primordial power spectrum is Gaussian, allowing us to simplify this 
expansion. This implies that 

Pi2 = 0, (132) 

due to 82 being symmetric under 5i — t- — At the next order, we find various contributions 
including the four-point function of 5i, which can be simplified using Wick's theorem. The final 
result for P13 and P22 will in general have two contributions: the contributions from IR modes 
and the contribution from UV modes, which we write in the following obvious notation 

Pi^{k) = PnMk, A) + Pn,uvik, A) (133) 
P22ik) = P22,iR{k, A) + P22,uv{k, A) (134) 

where we have reinstated the A dependence on the RHS as it separates the IR and the UV modes. 
By the IR contribution we mean the usual one-loop contribution, but cutoff at wavenumber A; 
this is given by the following loop integrals 

Pi3,/R(fc,A) = 3Pn(A:) y"^^Fg(q, -q, k)Pn(g) (135) 

P22,/R(fc,A) = 2 y""^ [Pi5(q,k-q)] Vn(g)Pn(|k-q|) (136) 

These contributions have a Feynman diagram representation that we present in Fig. O We have 
put a A superscript on the integrals as a reminder that they are cutoff by the smoothing function 
Wi^{k) that appears in Pn. By the UV contribution we mean the new fluid contribution that 
we are for the first time including in this work. This is given by the following 

{k,A) = F3,,{k,A)Pnik) 

(c2o(A) + c2o(A))P ^ ^ 
= - ' -Pii(^) (137) 

P22,UV {k,A) = AP22ik,A) (138) 

Here P13 is set by the (A dependent) sound speed and viscosity, and AP22 is set by the stochastic 
fluctuations that we elaborate on in Section [Dl the latter we find to be smaller than the former 
at low k as there is a suppression in the UV part of the integral by the transfer function. 
The IR contributions are associated with the long modes 61 running in the loop, while the UV 
contributions are associated with the short modes Sg running in the loop. 
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Figure 6: One loop Feynman diagrams for P13 (top) and P22 (bottom) after having extracted 
the time dependence. The crossed circles represent an insertion of the linear power spectrum 
and the loops represent convolution integrals cutoff at g ~ A. 

3.5 Cutoff Dependence of Fluid Parameters 

For the cutoff in the perturbative regime (A < /cnl); the A dependence of the fluid parameters 
is adequately described by the linear theory. This allows us to estimate the value of the fluid 
parameters and cf, and their time dependence as a function of the linear power spectrum. 
The sound speed is roughly given by the velocity dispersion, so in linear theory we estimate the 
sound speed by an integral over the velocity dispersion of the short modes. The linear theory is 
not applicable at very high k, so we shall include a constant correction as follows 



where is the velocity dispersion, a is an 0{1) constant of proportionality (which we will fix 
later), and Jj^dlnq = J (ilng(l — 1Va(^))^ since the sound speed arises from integrating out 
the short modes. In the A — )■ 00 limit, which we shall eventually take once we cancel the A 




(139) 
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dependence, we find that c^(A) is non-zero (due to the UV dependence), which we account for 
with the c^(oo) constant correction. 

Now we would like to relate the velocity to the power spectrum in terms of or 5. In the 
linear theory in a matter dominated universe universe, we have the growing mode solution 

Vi = -g^5L (140) 

which is the first order description of modes in the perturbative regime. This includes the short 
modes and 5s in the regime A < k < /cnl, which is a non-empty set if we choose small A. 

This gives the following linear relationship between the variance in the velocity and the 
density fluctuations 

where A|(/c) is related to the density power spectrum Piik) as given in eq. fl5^ . For a scale 
invariant primordial power spectrum we have 



Al{k) = ^A]{k) (141) 



k \2 



(see Fig. [3]) where we have taken into account the transfer function which separates the modes 
that enter before/after matter domination. Inserting this into eq. ( I139P leads to the following 
rough estimate for the sound speed 

where we have ignored logarithmic corrections, etc, so these estimates are only rough. In order to 
probe baryon-acoustic-oscillations, we shall need the effective description in the regime k > fcgg 
and so we need to take A ^ fceq which is the latter result. The full result from carrying out the 
integral over wavenumber is given in Fig. [TJ where we form a linear combination of pressure and 
viscosity. In the next subsection we will fix the coefficient a and the large A asymptotic value 
(c^(oo)) that were used to produce this plot. 

Notice that is time dependent, due to the l/'H? factor, and that it explicitly depends on 
the cutoff scale, due to the 1/ factor. A similar scaling goes through for the viscosity c^. This 
leads to the l/'H? oc a scaling that we stated earlier in eq. (11061) . 

3.5.1 Cutoff Independence of Physical Results 

The cutoff A explicitly alters the density field 5/, velocity field 9i, etc, and hence it affects the 
size of the loops. It also affects the size of the (bare) fluid parameters Cs,c^ etc, as summarized 
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Figure 7: The (bare) fluid parameter (c^(A) + /c^(A))/c^ at z = 0. The A dependence is chosen 
to cancel against the A dependence of the loop integral in P13. Note that the fluid parameter is 
non-zero as A — )■ 00, which we call c^(oo) + /c^(oo). The growth parameter / is different from 
1 in a ACDM universe which is described in Section I 



in Fig. [71 However, all the cutoff dependence must drop out of any physical results. In order to 
ensure this occurs, here we examine the form of the loop integrals f ll35|136p . 

By using isotropy, the full angular integral in Pizjr can be performed, and the azimuthal 
integral can be performed in P22 this leads to 
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2\2 



(144) 



(145) 



In order to extract the A dependence of Pi3(fc, A) we take the large r limit inside the integrand 
(which corresponds to the regime /c <^ A). In this limit the integrand approaches the value 



31 



-488/5, leading to 

Pi3Mk,A) = PisMk.M) - ^^l^^ii(^) tdqPniq) (146) 

where Ai is an arbitrary scale in the regime /c <^ Ai < A. On the other hand, the UV contribution 
is given in eq. ( I137P as 

(c2 (A) + c2 (A))P 

Notice that the A dependent parts of the IR and UV pieces have precisely the same time and 
/c-dependence ~ a^k'^Pii{k). Hence in order for the result to be explicitly cutoff independent, 
we must have 

<o(A) + <o(A) = (^^^^cigPn(g)) + <o(oo) + <o(oo) (148) 

This fixes the coefficient a that we mentioned earlier at a = ^gl^| (when combined with the 
viscosity term). The constant contributions are determined by explicit matching to numerical 
simulations. 

To make this more precise, let's separate the bare fluid parameters into a renormalized part 
^ren,o{''^ten) and a couuter-term c^^^o(^'-en, A) at some renormalization scale kren 

Clnfiikren) + cL,o(^ren, A) = ^(A) + o(A) (149) 

We define the counterterm to cancel the loop correction at the renormalization scale kren, i-e., 
the counterterm is defined through 



kren'^ctr,o{kren-, A) 



Pl3Mkren,A) ' ^ ^^^^^ ' PuiK-en) = (150) 



while the renormalized piece is defined such that the total power spectrum agrees with the full 
non-linear P{k) result at this renormalization scale, i.e.. 



P{Ken: t) = a\r)P,,{Ken) " 2a\r) ^^^^^^^if^ + a\T)P22{Ken) (151) 

y rLr 



Alternatively, by measuring the bare couplings directly from a measurement of the stress- 
tensor [t^-']k in simulations, we fix c^q(A) + /c^q(A) at some chosen A. The results of this 
numerical work we describe in detail in Section HI 
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3.6 Generalization to ACDM 



In the previous sections we have focussed on a matter dominated era, in which case the only time 
scale is set by Hubble. When we include dark energy this changes the background dynamics and 
the evolution. Assuming the dark energy is a cosmological constant, and operating at the linear 
level, the velocity field is related to the density fiuctuation by 

vl = dL (152) 

Here / is related to the growth function 

D(a)^lHSH(a)[jj^, (153) 

by 

/ . §^ (154) 

dm a 

This can be evaluated in terms of hypergeometric functions, which we do not reproduce here. 
So the corresponding relationship between the variances is 



= ^^"{k) (155) 



This leads to the fiuid parameters carrying the following time dependence 

A) = i--_c2 o(A), cl{a, A) = J-——cl^{K) (156) 

Jo ^0^0 J0TtQJ->Q 

(where the subscripts indicate the 2 = value, as before). 

Given the numerical solution for D (see the blue curve in Fig. (|8])) and the time dependence 
of the fiuid parameters, one can in principle construct the full time dependent non-linear solution 
perturbatively using Green's functions; however this is quite non-trivial as the Green's function 
is not known. To make proceed, we can make use of an approximation that is known to work 
reasonably well; we assume that the time dependence of the -n}^ order term is simply D{k,T)"' 
|14j . So to generalize the expansion fll02p for the matter dominated universe to the ACDM 
universe, we write 

oo 

5Kk,r) = J^DirTW (157) 

n=l 

oo 

eiik,T) = -nfY,D{rYW (158) 



n=l 
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The approximation is finalized by taking each 5n(k) to be the value in the D ^ a theory, i.e., the 
previously found solution for the matter dominated universe, which we denote with an "EdS" 
subscript. The corresponding approximation for the one-loop power spectrum is 

P(A;,r) = D{TfP,,{k) 

+ D{TY[2P,,{k) + P,,{k)\^^^ (159) 

where Pn{k) is the expression from eq. fllOSp . Piz{k) is the expression from fll35p . and P22{k) 
is the expression from (11361) . The reason this approximation works well is ultimately due to the 
fact that the dimensionless matter density VLm is approximately given by 

~ P (160) 

in a ACDM universe at all times. A full treatment of the time dependence was performed recently 
in [32] in terms of numerically evaluated Green's functions. These results can be compared and 
are found to be remarkably similar. The approximate results obtained provide useful analytical 
results and intuition, while the full Green's functions can be used for improved accuracy. 

As before, P13 includes a UV contribution from the fluid parameters. Generalizing the 
previous result from (I137p to the ACDM case, we have 

(c2 (A) + /oc2 (A))P 
Piw(A^,A) = - ' LJlnf -P^^(^) (161) 

which is of the form reported earlier in the introduction in eq. ([2]) (where we suppressed the 
detailed time dependence). The A dependence of the fluid parameters is given by 



c.,o(A) + /oc.,o(A) - 1— ^ 4^2 



^dgPn(g)) +<o(cx)) + /o<o(oo) (162) 



Evidently, the important combination is + fcf,, which is the value we measured and reported 
on in Fig. [71 



3.7 Summing the Linear Terms 

In the previous section we treated the fluid terms perturbatively. This meant we took the 
leading term in the expansion to be the usual growing mode in a matter dominated era 61 = 
f 61 oc D{a), and the fluid parameters provided corrections in an expansion in powers of the 
scale factor. However, since the sound speed and viscosity enter the linear theory, we can resum 
their contributions and form a new type of expansion in powers of the density field. In this way, 
all linear terms enter at first order, and only non-linear terms enter at second order, etc. This 
method treats the sounds speed and viscosity as independent parameters that can be measured 
separately. Although this is not necessarily possible in practice due to degeneracy with higher 
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order terms in the stress-tensor expansion, as we explain in Section I5.2[ this gives a sense of the 
consequences of summing a large number of terms. 
Let us write the expansion schematically as 

oo 

5i{k,T) = 5^5("Hk,r) (163) 

n=l 

oo 

9i{k,T) = 5^0(")(k,r) (164) 

n=l 

The equations of motion dnZEH]) give us the following linear order equations 

^ + ^«=0 (165) 
or 

ar 1 
_c2A;25(i) + i^^(i) = (166) 

By substituting eq. fll65p into eq. (I166p and re-arranging, we obtain a second order ODE for 
(5^^^. Lets express b^^^ in terms of a growth factor Dik^r) and a stochastic variable 5i(k), i.e., 

5W(k,r) = D(A;,r)5i(k) (167) 

The growth factor D satisfies the same ODE as namely 

We impose the asymptotic condition D{k,T) — ?■ a(r) for small a. The solution is plotted in 
Figure [8] for typical values of cj., c^. 

Although there is no simple analytical form for D{k,T), we can exhibit its structure. In 
particular, it has a self similar behavior, making it (up to a rescaling) only a function of a 
combination of a particular product of k, r, rather than k and r independently (the product is 
k T when the fluid parameters are treated as time independent and y/k t when they are treated 
as time dependent). For the case in which we take c^, cl to be time independent, let us rescale 
time to the following dimensionless variable 

T = ,/^,kT (169) 

We then find that in a matter dominated universe the ODE (11681) simplifies to 

d^D 1 ( dD 6 / 1T2 



dT"^ T \ A J dT T2 V b Q 

35 



b-r H^-Tp, l-T^ D = (170) 
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Figure 8: Re-summed growth function D{k,a) with k fixed at as a function of scale factor 
a, with standard cosmological parameters for a ACDM universe. The upper (red) curve is for 
k = 0, i.e, the usual growth function in a ACDM universe, which also coincides with SPT for all 
k with vanishing fiuid parameters. The lower (blue) curve is for the EFT with k = 0.2 [h/Mpc] 
and representative fiuid parameters = 7.2 x 10~^c^ and = 2.7 x 10~^c^ (although they 
exhibit degeneracy, as we discuss in Section [372]) . 

where b = Cy/cg. Of course this ODE has infinitely many solutions. Let us focus on one particular 
solution, which we denote T>{T), that satisfies the special asymptotic condition V{T) — )■ for 
small T. For any value of k the solution for D{k, r) is obtained from the one parameter function 
V{T) by 

gftr) = ^'^f> (171) 

This clearly has the correct asymptotic behavior D{k,T) — > a = (t/tq)^ for small a, since 
V{T) T2 for small T. 

In the case in which we take the fiuid parameters to be time dependent, as examined in 
Section 1331 with time dependence given in eq. fll06p . the analysis is slightly altered. In this case 
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we introduce the dimensionless variable 



y^(c,oC.,o)iv/fc^ (172) 



'To 

and the corresponding ODE in a matter dominated universe is 

where 6o = c^fl/csfi- Note the different powers of T in the parenthesis, compared to eq. fll7Up . 
Again we define the function D(T) as the solution to this ODE with asymptotic condition 
P(T) —7- for small T . The corresponding solution for the growth factor is 

( {csflC^fi)~i Vk t/ ^ ) 
= /c^fcr ^^^^^ 

For the power spectrum, we simply use the same form as before fll59p . but now dropping 
the Fa^i terms as they are built (and re-summed) into the linear piece. This approximation 
for P{k, t) is overly simplistic, however, since the growth function D is k dependent. A better 
approximation is to embed D{k,T) inside the convolution integrals of fll35|136p . Indeed we 
expect it to give somewhat accurate results, as has been the case in related calculations [H] and 
we shall use this approximation in Section HI The differential equations, whose solutions give 
the first few terms in the exact expansion, are provided in Appendix O 



4 Power Spectrum Results 

By matching to N-body simulations, as described in detail in Ref. [32], one can measure the 
linear combination c^q + /qC^^q. For simulation parameters: Qm = 0.25, Qa = 0.75, h = 0.7 
{H = 70 [km/s/Mpc]), dg = 0.8, and Ug = 1, with measurements described taking place at 
z = 0, and choosing a smoothing scale of A = 1/3 [h/Mpc] it was found 

<o + /o<o^9x 10-^c2. (175) 

Having obtained this linear combination, this completes the required quantities in order to 
compute the one-loop power spectrum that we derived earlier to the desired approximation. The 
single combination + fcl can be used as a single insertion by the formulae for P{k) derived 
in Section 13. 4[ or we can assume approximate individual values for cj. and separately in the 
re-summed formulae for the growth function D{k, r) and hence P{k) as indicated in Section IXTf 
although there is degeneracy in their values as we explain later in Section 15.21 

Our results for the full power spectrum P{k) in both of these approximations is given in Fig. |9l 
It is normalized to the no-wiggle power spectrum Pnw{k) of [37] is the hnear power spectrum 
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Figure 9: Power spectrum, normalized to the no- wiggle spectrum of [37] at 2; = for ACDM 
universe. Red is the linear theory (from CAMB), black is the full non-linear reference value, ma- 
genta is the one-loop SPT, green is the one-loop EFT with a single insertion of fluid parameters, 
and blue is the one-loop EFT with a type of re-summation of fluid parameters. 

without baryon- acoustic-oscillations. For convenience, we have taken the large A limit, by using 
the A dependence that we derived earlier in eq. (11621) . This causes the fluid parameters to 
asymptote to a slightly lower value, roughly 

4o(oo) + /o<o(oo) ^ 8 X 10-^ (176) 

In Fig. [H] we have also included the result for SPT for comparison and a non-linear reference 
value. We see that the power spectra of the EFT are much better than both the linear theory 
and SPT. The re-summed case is arguably better at lower /c, though it is somewhat degenerate 
with higher order effects. The single insertion is very accurate also. The fact that the fluid 
parameters asymptote to a finite non-zero value in the A — )■ cx) limit represents the finite error 
made in SPT. When made dimensionless, by multiplicity by c^k^ /H^ (say at 2; = 0), it leads to 
an important correction to the power spectrum as we increase /c, as seen in Fig. O Note that 
the EFT result is quite accurate; roughly at the 1% level out to ~ few x O.lh/Mpc. 
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5 Discussion 



In this section we briefly mention some interesting issues surrounding the effective fluid, including 
its Reynolds number, degeneracy of parameters, the inclusion of collisions or wave-like behavior, 
higher order moments, and the velocity field. 



5.1 Reynolds number 

For viscous fluids there is a famous dimensionless number which captures its tendency for laminar 
or turbulent flow; the "Reynolds number" . The Reynolds number is defined as 

fl. . (177) 
r] 

where rj is shear viscosity, p is density, w is a characteristic velocity, and L is a characteristic 
length scale. This is 

Hence the Reynolds number is not very large, and the system is therefore not turbulent. Fur- 
thermore, if we were to estimate the viscosity by Hubble friction, then we would have Re ^ 5 
and so the Reynolds number would be even smaller in the linear or weakly non-linear regime. 

For cosmological parameters pf, ~ 3 x 10"^*^ [g/cm^], H = 70 [km/s/Mpc], and if we take a 
plausible value for the shear viscosity of c^^ ~ 2 x 10~^c^, then the viscosity coefficient is found 
to be 

r/~20Pas (179) 

which is perhaps surprisingly not too far from unity in SI units. (For instance, it is somewhat 
similar to the viscosity of some everyday items, such as chocolate syrup.) A proper measurement 
of cj.^ would come from a detailed measurement of vorticity; a point that we will return to in 
the following subsection. 



5.2 Degeneracy in Parameters 

Earlier we discussed the individual parameters: the sound speed cj., shear viscosity c^^, and bulk 
viscosity c^„. We showed that by taking the curl of the Euler equation we obtained an equation 
for vorticity that involves Cg„. Hence a careful analysis of vorticity could reveal the value of 
c^^ - such a value would enter our discussion of "Reynolds number" of the previous subsection, 
although not the bulk of this paper. Since the vorticity is rather small, this would be non-trivial 
to measure, though possible. 

On the other hand, by taking the divergence of the Euler equation we obtained coupled 
equations for 6i and 6i that involves the sound speed and the combination of viscosity = 



39 



+ c^j,. One could try to measure these two parameters independently using eqs. (pOllGTl) 
from N-body simulations. However one should be careful as to how to interpret this result. At 
the one-loop level, we saw that it was only a certain linear combination that appeared in the 
result, namely + / c^. In other words, the two parameters appear in a degenerate way at one- 
loop. This degeneracy would be broken at higher loop order. However, to be self-consistent one 
should then also include new couplings (for instance, representing higher derivative operators in 
the stress-tensor expansion) which would also enter, leading to a new constraint and a new type 
of degeneracy with the new parameters. 

This reason for this degeneracy is the following: In a universe in which one could track the 
full evolution of the initial state, one would observe that and affect the one-loop theory 
differently. However, in the real universe, there is a growing mode and a decaying mode. In 
practice, one does not track the decaying mode, only the growing mode, as studied here in this 
paper. For this single mode the parameters enter in a special linear combination. 



5.3 Interactions 

In this paper, we have treated dark matter as being comprised of collisionless particles, interact- 
ing only through gravity. Of course we expect that there are also some finite non-gravitational 
interactions. These effects can be treated perturbatively in the effective field theory by identi- 
fying the relevant length scale, which is the mean free path between scattering. 

As an example, for WIMPs the scattering cross section is roughly a ~ g'^/m'^r, where (7 ~ 0.1 
is a dimensionless coupling and mw ~ 100 GeV is of order the weak scale. The mean free path 
between scatterings is 

Ampp = — ~ ^ (180) 
The background density, in Planck units, is pb ~ 10~^^°Mpj. Hence, 



A. 



io"^Hf-Vr y (181) 

V ^ / VlOOGeV/ ^ ' 



where dn = ^/Hq is the present Hubble length. As another example, for QCD-axions scattering 
due the (p^ term in the potential V{(j)) = A|^^(l — cos{(f) / Fpq)) , the mean free path is much 
larger still. 

On the other hand, gravity introduces a non-linear scale of the order A^l ~ 10"^ cIh (which 
can be thought of as the mean free path between gravitational scattering of a point particle 
off a large non-linear clump). Since the collisional mean free path in these examples satisfies 
Amfp ^ Anl it can be ignored at first approximation. Though in principle it can be included 
perturbatively in the effective field theory, but suppressed by a hierarchy Anl/Amfp, which may 
be of interest for some highly non-standard dark matter candidates. 
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5.4 Wave- like Behavior 



In this paper we have treated the dark matter as comprised of classical point-like particles. This 
obviously ignores its quantum mechanical wave-like behavior. For most dark matter candidates, 
such as a typical WIMP with a weak scale mass, the de Broglie wavelength is extremely small 
and ignorable. For extremely light (pseudo)-scalars, such as axions, it is conceivable that their 
de Broglie wavelength is large and relevant. 

For a classical scalar field and also for a Bose-Einstein condensate [31], one can show that in 
the linear theory, there is a correction to the pressure of the form 

= -J^^'P (182) 

this provide a contribution to a type of scale dependent sound speed ^c^ ~ Tv^k"^ / {a?"m?). Now 
recall that the characteristic correction from the sound speed is ~ c^^^/ {H'^a?). This means that 
a rough estimate for the dimensionless correction from the quantum character of the particles is 
(at ^ = 0) 

quantum correction ~ — - — ^ (183) 

The relative size of this contribution obviously depends on the mass of the particle m. It is 
important to note that in the point-particle treatment, the mass m dropped out of all results. 
But by including such UV physics, we gain more sensitivity in the effective field theory to such 
physical parameters. 

The dark matter particle mass m can in principle be very small. For instance, in the so-called 
string axiverse it is suggested that there may be a range of extremely light axions [SB]; one of 
which could provide the bulk of the dark matter (though there are important constraints from 
isocurvature bounds on light axions [391 HQ], while the classic axion window is still very promising 
[m Il2])- Here we would like to mention that the mass of a dark matter particle presumably 
cannot be arbitrarily small because its de Broglie wavelength A^b ~ hj {mv) would then smear 
it out over scales larger than that of a galaxy Lgai, and yet we know dark matter clumps on 
galactic scales. By imposing XdB < Lgai this gives the bound 

m > (184) 

LgalV 

Hence we have a bound on the dimensionless correction from the wave-like character of light 
scalars as 

quantum correction < — (185) 

where w is a characteristic dispersion velocity associated with the dark matter. By estimating 
V ~ 10~^c, then we can estimate Hq/v ~ k^L-, leading to the rough bound 

quantum correction < ^ (186) 
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Since Lgai is much smaller than the non-linear scale (for instance, Lgai may be as small as dwarf 
galaxy size) we see that the quantum correction must always be very small in the regime in 
which the effective field theory is valid (i.e., k < k^L). 

5.5 Higher Order Moments 

In principle, one can study higher order moment of the Boltzmann equation. In Section [2] we 
considered the zeroth moment (continuity) and first moment (Euler), and then built a derivative 
expansion for the effective stress-tensor that appears on the right hand side of the Euler equation. 
The stress-tensor involves two contributions: kinetic and gravitational. The kinetic piece is 
includes the second moment of the velocity distribution (minus the long modes), and so it evolves 
under the second moment of the Boltzmann equation. The trace of the kinetic part of the stress- 
tensor is proportional to a type of "kinetic temperature" T ~ ^i/pi- Although the system is not 
in thermal equilibrium, so this name is only by analogy to classical systems which are. 

However, since the stress-tensor also includes the gravitational piece we do not have an 
evolution equation for the full stress-tensor. Parametrically these two contributions are of the 
same order (for instance they cancel each other in the virial limit). So this requires the use of 
a derivative expansion to capture the effects of these higher order moments and interactions, in 
the effective field theory sense. 

5.6 Velocity Field 

Let us also make some comments on the computation of correlation functions, such as the two- 
point, involving the smoothed velocity field. In Section [3.5.11 we demonstrated how the cutoff 
dependence in the expansion for Si cancels out when we form the two-point correlation function 
for density, which involved the fluid parameter's canceling the A dependence of the P13 loop. 
However, once the fluid counter-terms are introduced to cancel this dependence, then they cannot 
also be used to cancel the cutoff dependence appearing in the loops of correlation functions of 
other fields. In particular, consider the velocity field. Recall that it is defined by the ratio 
vj = T^ll pi- If we Fourier transform this, then examine the regime /c < A, we are still left with 
A-dependence, even though both pi and tt^ are cutoff independent for k < A. This means that 
the velocity field vj is inherently cutoff dependent even in the k <^ A regime, while 61 is not. 

6 Summary and Outlook 

In this paper we have examined and developed the effective field theory of dark matter and 
structure formation on sub-horizon scales, emphasizing detailed analytical and semi-analytical 
results, including a recursion relation for the perturbative expansion and an approximate ex- 
traction of the time dependence of the growing modes. This works compliments and extends 
the important recent work in Ref. [32]. These works can be viewed as a precise realization of 
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the conceptual foundation laid out in Ref. [31], where special focus was placed on the issue of 
back-reaction at the scale of the horizon, though the present focus is on sub-horizon scales and 
the explicit computation of the power spectrum. The effective field theory is an expansion for 
wave numbers k less than the cosmological non- linear scale fc^L. It is a cosmological fluid descrip- 
tion for cold dark matter, and by extension all matter including baryons which trace the dark 
matter. The microphysical description was in terms of a classical gas of point particles, which 
we smoothed at the level of the Boltzmann equation and used the Newtonian approximation for 
sub-horizon modes. We exhibited the various couplings that appear in the effective field theory, 
namely pressure and viscosity, whose linear combination was obtained by matching to N-body 
simulations in Ref. [32], finding cl + f cl ~ 10~^c^. This represents the finite error made in 
standard perturbation theory and has important consequences for the power spectrum. We see 
that standard perturbation theory would only be correct if the linear power spectrum was very 
UV soft so that the theory remains perturbative to arbitrarily high k; this would allow one to 
send A — oo and there would be no stress-tensor at all as it is sourced only be the short modes. 
However the presence of the non-linear scale forces one to introduce the cutoff A and a finite 
stress-tensor, which evidently does not vanish in the large A limit (which can be formally taken 
order by order after the A dependence is cancelled). 

We developed the perturbative expansion for the power spectrum, which we recast into a 
recursive formula and extracted the time dependence in a convenient way. The power spectrum 
was then computed at the one-loop order. We found that the corrections from the fiuid param- 
eters led to a power spectrum in good agreement with the full non-linear spectrum. Unlike the 
standard perturbation theory that deviates substantially from the true non-linear power spec- 
trum, especially at low z, suggesting that the non-linear wavenumber is low, the effective field 
theory exhibits ~ 1% level accuracy for k ~ few x 0.1 h/Mpc, suggesting that the true non-linear 
wavenumber may be higher than ordinarily thought. Furthermore, the general success of the 
effective field theory approach suggests that any deviations from the ACDM model should fit 
into this framework by altering the fiuid parameters. 

The effective field theory approach to large scale structure formation is complimentary to 
N-body simulations by providing an elegant fiuid description. This provides intuition for various 
non-linear effects, as well as providing computational efficiency, since the numerics required to 
measure the fluid parameters can be less computationally expensive than a full scale simulation. 
Of course, since the couplings are UV sensitive, it still requires the use of some form of N-body 
simulation to flx the physical parameters, either by matching to the stress-tensor directly or 
to observables; a point analyzed in detail in Ref. [32]. But this matching is only for a small 
number of physical parameters at some scale and then the constructed fleld theory is predictive 
at other scales. The formulae for the power spectrum, exhibited in eqs. fll31H137tll48p for matter 
dominated and eqs. f ll59H161^ for ACDM, should be quite powerful and convenient in this area 
of cosmological research. 

There are several possible extensions of this work. A flrst extension is to go beyond the 
one-loop order to two-loop, or higher. It is important to note that the treatment can in principle 
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capture the full power spectrum to arbitrary accuracy if carried out to the required order. Here 
we have computed the power spectrum at one-loop, which is order S^, but higher order is possible. 
This will require the measurement of several new parameters that will enter the effective stress- 
tensor at higher order. Another extension is to fully measure the stochastic fluctuations, which 
will involve measuring the correlation function of the stress-tensor with itself. These effects 
are somewhat reduced at low k, especially due to the suppression of modes in the integrand 
due to the turnover in the transfer function in eq. (11381) . but are of significant interest and 
should improve the agreement at higher k. Another extension is to include the small but finite 
contributions from vorticity, or to compute higher order N-point functions, which can probe 
non-Gaussianity, or to consider different cosmologies, and to include baryons, etc. 

In general it is essential to gain insight and precision into the mapping between the micro- 
physics that determines the early universe, including the distribution of primordial fluctuations 
and the contents of the universe, and the output universe that we can observe today. This 
approach, when complimented with N-body simulations and observations, may ultimately give 
new insights into fundamental questions in cosmology. This is an exciting and promising way to 
learn about fundamental physics. 
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Although we use the full stress-tensor in fl33l) - fl35|) as a generating functional of the effective 
theory, in this appendix we demonstrate that we can separate out the long modes from the short 
modes in the Euler equation. To do so, we define the short modes to be 



A Short Modes 




m ^ /o?3p(p'-pi(x))(p^' -p/(x))/(x,p) 
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where p\ (x) = mav] (x) . Note that a^^ 7^ a^^ — a]'' , but they are related as follows 

- + M^/Ia + - PrnVl) + Vl{7^^ - PrnV^)] ^ (192) 



The second term is approximately pivlvj (so it approximately cancels with —pivlvj in nf) and 
the final term is small (as it is an overlap between short and long modes). Following the methods 
of [31] we obtain 



1 

A4 



(193) 



Similarly, one can prove that satisfies 



So altogether we obtain the effective stress-tensor 
where 



A 87rGa2 
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(196) 
(197) 



We see that [t^-'Ja is sourced by short wavelength fluctuations plus higher derivative corrections. 
Note that by taking the derivative dj this leading piece becomes 



5.Ma = 5.K']a + M 



with 



[psdi(t)s]k = ^ma ^di(j) sfiMWAix - x„) - [pid., 
where the flrst term in fll99p is given by 
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1^a(x„-x,0 iyA(x-x„) (200) 
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and the second term in fll99p can be expanded as 

[pidi(t)s]K = -^Pidid^(pi + ... (201) 

and this term should be included since it involves the background piece pb, and so it includes a 
first order contribution. 

B Trace of Stress- Tensor 

It is of some interest to compute the trace of the stress-tensor, which is the so-called "mechanical 
pressure" . This includes the gravitational piece 

». = -^ + |^ (202) 
The first term is approximately given by 

^ I d'^' iyA(x - x') [5p(x')</.(x') - ^ 5p„(x')</'n(x')] (203) 

1 y J^"^' -M|x„-x.|^ _ ) 



87rGa2 2 



where we have used the identity 



(V0)2 = -0V^0 + ^V2(02) (205) 



and dropped all terms that are suppressed by the ratio of low A;- modes to high /c- modes in f l203p . 
If we let an s-subscript denote the short modes (see Appendix [A] for details), then the trace of 
the stress-tensor is roughly 

[r]A ^ f rf^x' W^{^ - x') [p(x')(t;.(x')2 + ^0s(x )) Ps,n(x' )0n(x )] (206) 

The background pressure has the zero mode contribution 

Vi> = \\rV) (207) 

where we have ignored a correction from the bulk viscosity. There are also stochastic fluctuations 
to the pressure, which we discuss separately in Appendix [Dl Now, since the density field p(x) can 
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be arbitrarily large for dense objects on small scales, it suggests that each of the contributions to 
the renormalized pressure, both the kinetic and the gravitational, can be quite large. However, 
for virialized structures, these two terms cancel each other [31]. Hence the only significant 
contribution to the integral comes from modes of order k ~ /c^l which have yet to virialize. For 
a pure Einstein de Sitter universe, this leads to the estimate pb ~ 10~^pi,c^, where the factor of 
10~^c^ is the typical value of the potential from the primordial power spectrum. However, due 
to the turn-over in the power spectrum at matter-radiation equality, we are led to an value that 
is about an order of magnitude smaller; see Section 13.51 



C Exact Expansion 



The expansion in Section 13.71 should be compared to the exact expansion, which we briefly 
mention. The second order density fluctuation 5*^^^ can be written in terms of a pair of time 
dependent kernels D2A and D2B as follows 



5(2) (k,r) 



X 



|^5,(k-k05,(k') 

5a(k, k')Z^2A(k, k', r) + 2/3(k, k')Z?2s(k, k', r) 



The kernels satisfy the following ODEs 

LD2A 
LD2B 
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5 

7dD 



^n'(n^-'^)D{\k-k'\)D{k') 



3^2 



(Ik-k'l)^(A;') 
2 dr dr 



(208) 

(209) 
(210) 



where L is the linear operator 



L = 



dT^ 



U 1 



d 3^ ,9 



2c^2 



(211) 



We have suppressed the r dependence in the argument of D and the r, k, k' dependence in 
the argument of D2* in f l209ll210|) . As usual we impose the asymptotic condition D2* — ^ 0? for 
small a. By numerically solving this pair of ODEs, we can compare to the approximation in 
fll581ll59| in which we replace D2*(k, k',r) — t- D{k,T)'^ in the power spectrum. At third order, 
we again expand 6^^^ in kernels Ds* to obtain ODEs whose solutions can be compared to the 
approximation in f ll58lll59| in which we replace D3^,(k, k', k", r) — > D{k, r)^. These expressions 
can also be given in terms of Green's functions as we emphasized recently in [32]. As can be 
checked, the approximate and intuitive analytical results are quite close to the exact results. 
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D Stochastic Fluctuations and Pressure 



The effective stress-tensor [r*-']A ffuctuates with the short modes, leading to a stochastic depar- 
ture from the results computed thus far. This effect should be reduced at low k since the transfer 
function softens the UV modes in the corresponding loop integral, but we would like to mention 
the formal procedure to include such finite size effects here, although we will not compute it 
precisely. 

Again ignoring vorticity, but allowing for stochastic fluctuations, the evolution equations for 
the pair 6i and 9i are as given in eqs. (PTIIUS]) . but with an additional stochastic source term 

^ + (^1 = -y^«(k,k')5Kk-k')^Kk') (212) 

^ + m + ^n'Q^Si = - / 7^/3(k, k%{k - k')9iik') + cy6i - - J, (213) 

ar Z J [Zir)-^ H 

where js = js(k, r). In position space js is defined from a scalar contraction of the fluctuations 
in the stress-tensor, namely 

js = -d,d,AT'^ (214) 

Pb 

where Ar*-^ is implicitly defined through eq. ( |4T|) . So (up to higher order corrections) it is 
therefore related to the function Ai by 

a'A = cld'6,~'^+js (215) 

Note that js does not enter the two-point correlation functions we defined earlier, such as {Ai 6i), 
since the averaging annihilates any overlap between the non-stochastic pieces Si, 9i and the 
stochastic piece js- However, a two-point correlation function involving only the stochastic piece 
{jsjs) would be non-zero. 

We treat the js term as entering at higher order in a field expansion of the form fll63lll64| . As 
before, the first order term 5^^^ is given by the growth function D, defined through eqs. (1167111681) . 
At second order, let us split the density field into two pieces 

5(2) = 5(2) + ^^(2) ^216) 

where 5q is defined as the contribution that arises even in the js — limit that we computed 
previously, i.e., is given by eq. f l208p (or earlier approximate forms). While A^''^) jg the 
correction that arises from the stochastic contribution js. This new piece satisfies 

LA(5(2)=j, (217) 
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Formally, this may be solved by Green's functions. So let us introduce G{t), defined as a 
particular solution of the same differential equation, but with a Dirac-delta function RHS 

LG{k,T,T') = 6Ur -r') (218) 

Then the contribution to the second order density from the stochastic fluctuations may be 
formally given by 

A5(^)(k,r) = Jdr'G{k,r,r')Uk,r') (219) 

We note that one can use this Green's function to express the solution in f l208p . 

The two-point correlation function {5^'^^ 5^^^) has three contributions when we expand using 
fl2T6|) . The first term (^J^^^^^) 

is as we computed earlier. The cross term ((5o^'' A^*-^-*) vanishes, 
leaving only a single new term (A^*-^-* A^^^^), i.e., 

(5(2) 5(2)) = (^(2) ^(2)) + (^^(2) ^^(2)) (220) 

The new term may be expressed in terms of the Green's functions and the stochastic contraction 
of the stress-tensor js as 

(A5(2)(k,r) A5(2)(k',r)) 

= I dr'l dr" G{k, r, r') G{k', r, r") (j,(k, r')j.(k', r")) (221) 

This requires one to obtain the ensemble average (js(k, r')js(k', r")), which may be difficult to 
obtain numerically. Using isotropy we may partially simplify it to 

(j,(k, r')j.(k', t")) = {2nf6U^ + k')P,,{k, r', r") (222) 

This leads to the following contribution to the two-point power spectrum 

AP22{k,T) = Jdr' JdT"G{k,T,T')G{k,T,T")Pjj{k,r',r") (223) 

This term should be understood as providing the UV part {k > A) of the P22 Feynman 
diagram drawn in Fig. [61 We have checked that these corrections are small at low k due to the 
transfer function T{k) in the real universe (though this would not be the case in pure Einstein de 
Sitter), and will therefore was not included in our numerics. On the other hand, the contribution 
from the UV part of P13 contains the leading departure from SPT and is connected to the fluid 
parameters, as we have computed. 

Related to this is how much the effective stress-tensor [t^-'Ia varies from patch to patch. This 
is measured by contractions of the variance 

Var[r]f ^-'(t, t') = {[r%{t) [T^'%{t')) - {[T%{t)) {[T^'%{t')) (224) 
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If we contract over i,j and and set t' = t, then this measures the fluctuations in the 

pressure; a type of stochastic pressure 

i^pf = I mi) - (Ma)^) (225) 

This is the fluctuations around the background value given by the mean pb = |([t]a) fl207p . The 
absolute pressure does not affect the Newtonian dynamics (although it does affect the metric in 
GR). Instead the important fluctuations are those that arise from derivatives of the stress-tensor 
and are connected to the UV part of the P22{k) contribution to the power spectrum. These are 
small at low k, due to the transfer function in high /c-modes as we mentioned previously, but 
can play an important role at higher k approaching kisf^. 
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